Failure of random matrix theory to correctly describe quantum dynamics 
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Consider a classically chaotic system which is described 
by a Hamiltonian Tlo- At t = the Hamiltonian under- 
goes a sudden-change Ho i-^ H. We consider the quantum- 
mechanical spreading of the evolving energy distribution, and 
argue that it cannot be analyzed using a random-matrix the- 
ory (RMT) approach. RMT can be trusted only to the extend 
that it gives trivial results that are implied by first-order per- 
turbation theory. Non-perturbative effects are sensitive to the 
underlying classical dynamics, and therefore the Ti — > be- 
havior for effective RMT models is strikingly different from 
the correct semiclassical limit. 



Consider a system whose total Hamiltonian is 7i = 
H{Q, P; x), where {Q,P) is a set of canonical coordi- 
nates, and X is a constant parameter. Wc assume that the 
preparation and the representation of the system are de- 
termined by the Hamiltonian Hq — 'H{Q, P; xq), and that 
both Ho and Tl generate classically chaotic dynamics of 
similar nature. Moreover, we assume that dx = (x—xq) 
is classically small, meaning that it is possible to apply 
linear analysis in order to describe how the energy sur- 
faces H{Q, P;x) — E are deformed as a result of changing 
the value of x. Physically, going from 7-^o to TL may sig- 
nify a change of an external field, or switching on a per- 
turbation, or sudden-change of effective-interaction (as 
in molecular dynamics). Quantum mechanically, we can 
use a basis where = Eq has a diagonal representation, 
while 
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For reasonably small it follows from general semiclassi- 
cal considerations Q , that B is a banded matrix. Gener- 
ically, this matrix looks random, as if its off-diagonal el- 
ements were independent random numbers. 

It was the idea of Wigner |^ forty years ago, to study 
a simplified model, where the Hamiltonian is given by 
Eq. (|l|), and where B is a Banded Random Matrix (BRM) 

. This approach is attractive both analytically and nu- 
merically. Analytical calculations are greatly simplified 
by the assumption that the off-diagonal terms can be 
treated as independent random numbers. Also from nu- 
merical point of view it is quite a tough task to calculate 
the true matrix elements of the B matrix. It requires a 
preliminary step where the chaotic TYq is diagonalized. 
Due to memory limitations one ends up with quite small 
matrices. We can think of Eq. (P as describing ficti- 
tious motion on a lattice. For the model below (Eq. ||) 
we were able to handle iV = 5000 sites maximum. This 
should be contrasted with BRM simulations, where us- 



ing self-expanding algorithm [Q we were able to handle 
N = 100000 sites along with significantly reduced CPU 
time. 

However, the applicability of the RMT approach is a 
matter of conjecture. Obviously this conjecture should 
be tested. To be more specific, one should be aware that 
there is a hierarchy of challenges where the applicability 
of the RMT conjecture should be tested. Namely: The 
study of spectral statistics; The study of eigenstates; The 
study of quantum dynamics. While the issue of spec- 
tral statistics has become a major subject in "quantum 
chaos" studies |^ , the two other issues are barely treated. 
In a previous study Q we have demonstrated that the 
RMT approach is capable of giving the right qualitative 
picture of the parametric evolution of the eigenstates. As 
5x is increased the eigenstates of Eq. (P change in a qual- 
itative agreement with Wigner 's theory. Still, RMT fails 
to capture non-universal system-specific features. 

In this Letter we turn to the study of quantum dy- 
namics. Here we are going to end up with a much more 
alarming claim. Namely, the RMT approach fails to give 
the correct dynamical picture. RMT can be trusted only 
to the extend that it gives trivial results that are implied 
by first-order perturbation theory. Non-perturbative ef- 
fects are sensitive to the underlying classical dynamics, 
and therefore the ?i — > behavior for effective RMT mod- 
els is strikingly different from the correct semiclassical 
limit. In this Letter we are going to establish the failure 
of RMT for the case of dynamics which (for t > 0) is gen- 
erated by a time independent Hamiltonian. This we hope 
paves the way towards making an analogous statement 
regarding the response of driven systems ^ . 

In order to test the RMT conjecture, we are going to 
use the same 'direct' approach and the same model as 
in Ref. 0. For the first time we are going to compare 
the dynamics which is generated by a 'physical' Hamilto- 
nian, with the corresponding dynamics that is obtained 
from an effective BRM model (EBRM). The latter is con- 
structed by taking the matrix B of the 'physical' Hamil- 
tonian, and then randomizing the signs of its off-diagonal 
elements. Such operation destroys any correlations be- 
tween the matrix elements of B, while keeping the band 
profile unaffected ||^. We study the Hamiltonian 



H{Q,P;x) ^ hjPf+P^ + Ql+QD+x-QlQ 



(2) 



with x = Xq + Sx and xq = I. This Hamiltonian de- 
scribes the motion of a particle in a 2D well (2DW). The 
units are chosen such that the mass is equal to one, the 
frequency for small oscillations is one, and for Sx = 
the coefficient of the anharmonic term is also one. The 
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FIG. 1. (a) The classical energy spreading as a function of 
time, (b) The QM spreading for the 2DW Hamiltonian; (c) 
The QM spreading for the EBRM Hamiltonian. The energy 
in these simulations is _E ~ 3, and Sx = 0.2123. In (a) we see 
a crossover from ballistic spreading {SE oc t) to saturation 
{6E ~ const). Only one time scale {t^i ~ 1) is involved. The 
light dashed line has slope 1 and is drawn to guide the eye. In 
(b) we see that the classical behavior is approached as ?i ^ 0. 
In (c) we see the opposite trend: as ft ^ an intermediate 
stage of diffusion {SE oc y/t) develops. (Here the light dashed 
line has slope 1/2). Different lines correspond to different 
values of h as in (b), and additional curves {h = 0.009, 0.005) 
have been added. 

energy E is the only dimensionless parameter of the clas- 
sical motion. Our numerical study is focused on an en- 
ergy window around E ^ S where the motion is mainly 
chaotic. Upon quantization we have a second dimension- 
less parameter, which is the scaled h. Associated with 
h are two energy scales. One is the mean level spacing 
A (X with d = 2, and the other is the bandwidth 
Aft oc h. The second scale is further discussed below. 

It is useful to define a fluctuating quantity J-{t) = 
— {dH/dx) which for the 2DW model equals !F{t) ~ 
—Qi{t)Q2{t)- The auto-correlation function of !F{t) is 
denoted by C(r). The associated correlation time is de- 
noted by Tel. The power spectrum of the fluctuations 
C{uj) is the Fourier transform of C(t). The band pro- 
file of the matrix B satisfies the semiclassical relation 
|B„„|2 « {A/i27rh))C{{E,,~E^)/h). See Fig.2 of § for 
numerical demonstration. It is implied by this relation 
that the bandwidth is — 27r?i/r^i. For the quantum 
mechanical simulation the exact matrix B has been cal- 
culated numerically. Memory constraints limit the maxi- 
mum size (N) of the matrix that we can get. The EBRM 
Hamiltonian is obtained by randomizing the signs of the 
off-diagonal elements. A second, more 'loose' strategy, is 
to generate the EBRM B from scratch using the semi- 
classical band-profile as an input. The advantage of the 
latter strategy is that it opens the way to EBRM-model 
simulations with smaller h, where the required N is much 




FIG. 2. The width A x N{t) as a function of time. Differ- 
ent lines correspond to different ?i-values as in Fig. lb. The 
heavy dashed line is the classical 5E{t). Having separation of 
scales {AxN{t) ^ 5E{t)) is an indication for having a per- 
turbative spreading profile. The upper panel is for the 2DW 
Hamiltonian while the lower panel is for the EBRM model. 

larger. We have verified that the latter strategy gives 
numerical results that agree with the sign-randomization 
approach. 

The initial preparation is assumed to be microcanoni- 
cal. This means, in the classical case, an ergodic distribu- 
tion of initial 'points' on the energy surface TLq{Q, P) — 
E with ~ 3. In the quantum mechanical case we 
start each simulation with an initial eigenstate m that 
has an energy Em ~ E where 2.75 < E < 3.2. The 
time dependent evolution is determined by Schrodinger 
Equation. The probability distribution after time t is 
Pt{n\m). An average over initial state (m) is taken in or- 
der to get the average profile Pt {n — m). We characterize 
the evolving distribution using three different measures. 
The variance is M{t) = J2r ''^'^ ^t(r) , or in energy units 
it is 6E{t) = A X yjM{t). The width N{t) is defined as 
the r region that contains 50% of the probability. In case 
that we have a spreading profile that is characterized by 
a single energy scale, it is implied that N{t) and y/M{t) 
would be the same (up to a numerical factor). The sur- 
vival probability is P{t) — P((r=0). The results of the 
simulations are presented in the Figs 1-4. The analysis 
of these results is discussed below. 

Taking 7i to be a generator for the classical dynamics, 
the energy E{t) = Ho((5(i), P{t)) fluctuates. The fluctu- 
ations are characterized by the correlation time t^,, and 
by an amplitude SE^i . The initial preparation is assumed 
to be a microcanonical distribution that is supported by 
the energy-surface HoiQ^P) = E{0). For t > 0, the 
phase-space distribution spreads away from the initial 
surface. 'Points' of the evolving distribution move upon 
the energy-surfaces of Ti.{Q, P). We are interested in the 
distribution of the energy E{t) of the evolving 'points'. 
It is easily argued that for short times this distribution 
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FIG. 3. The width N{t) and the survival probabihty P{t) 
for the 2DW modeL Different Unes correspond to different 
R-values as in Fig. lb. Having N{t) = 1 or P{t) ~ 1 is 
the indication for having a standard perturbative spreading 
profile. 

evolves in a ballistic fashion. Then, for t ^ t^, due to 
ergodicity, a 'steady-state distribution' appears, where 
the evolving 'points' occupies an 'energy shell' in phase- 
space. The thickness of this energy shell equals 6E^i. 
Thus we have a crossover from ballistic energy spread- 
ing to saturation. The dynamics in the classical limit 
is fully characterized by the two classical parameters t^i 
and SE^,. 

A quantitative description of the classical spreading is 
easily obtained. A straightforward derivation leads to the 
following result for the spreading 

SE,,{t) ^5xx v/2(C(0) - C{t)) (3) 

As a particular result we get SE^i = SE^i{oo) ~ 
(5x^2C(0). The calculation of SE^i{t) for the model 
Hamiltonian is presented in Fig. la. It is implied by 
Eq.(||) that the spreading SE{t), from semiclassical point 
of view, is just a property of the band-profile. Thus, 
one may get to the wrong conclusion that models with 
the same band-profile should lead to the same SE{t), pro- 
vided the off-diagonal elements look random. If this were 
the case, it would be implied that the EBRM model 
would be equivalent to the 2DW model as far as the 
spreading SE{t) is concerned. Looking on Fig.l we see 
that this is not the case. As ?i ^ the EBRM model fur- 
ther and further deviates from the (correct) semiclassical 
expectation. 

In order to understand the observed results we would 
like to recall some of the theory of We already 

said that upon quantization we have the two energy 
scales A a h'^ and and A;, c>c 7i. Actually there is 
also a semiclassical energy scale Age c>c h'^^^. Associ- 
ated with these energy scales are three parametric scales 
Sx"^"" <^ feprt ^ Sxsc, where the strong inequalities 



hold in the h ^ limit. For the 2DW model, assum- 
ing ~ 3 we have the estimates 6x1"' ~ 3.8 * h^^'^ 
and 6xprt ~ 5.3 * h and 5xsc ~ 4 * Ti^^"^. In the stan- 
dard perturbative regime {Sx < 6x1"^) the eigenstates 
of Eq.(|l|) have a simple perturbative structure. In the 
extended perturbative regime (fe^" < 6x < Sx-p^t) the 
eigenstates of Eq. (||) have a core-tail structure that can 
be regarded as a generalization of Wigner's Lorentzian. 
In the non-perturbative regime {Sx > Sx^^t) the eigen- 
states have a purely non-perturbative structure. Depend- 
ing on whether it is the EBRM Hamiltonian or the 2DW 
Hamiltonian, this 'ergodic' non-perturbative structure 
is either semicircle-like or semiclassical-like respectively. 
The semiclassical regime (Sx > Sxsc) is contained in the 
non-perturbative regime. It is only there that we can 
trust detailed quantal-classical correspondence (QCC). 

For the purpose of the present analysis it is con- 
venient to specify the different regimes by regarding 
?i as a free parameter. Thus, the standard pertur- 
bative regime is ?i > Ccq„, the extended perturbative 
regime is Cp^ < h < Ccq„, and the non-perturbative 
regime is ?i < Cp,.t. The latter contains the semiclassi- 
cal regime Ti < Csc- Thus the semiclassical limit fi — > 
is a non-pcrturbative limit. In case of the 2DW 
model, the classical quantities are Ceq„, = 0.41 * dx^^^ 
and Cprt = 0.19 * dx and Csc = 0.12 * dx^/^. We have 
used in most of our numerical simulations Sx ^ 0.2. 
Larger Sx may take us out of the classical linear regime. 
For this value of Sx we get C<.q„ = 0.14 and Cprt = 0.04 
and Csc = 0.01. The smallest h value that we could allow 
without having memory-overflow was h = 0.015. This 
means that we were able to access the non-perturbative 
regime, though the semiclassical regime was out of reach. 

As explained in [|lO| the essential features of the spread- 
ing behavior in the perturbative regimes can be analyzed 
using first order perturbation theory (FOPT). Since cor- 
relations between off-diagonal elements are not impor- 
tant for FOPT, it follows that the EBRM-model and 
the 2DW model should be trivially equivalent in such 
case. It is only in the non-perturbative regime where 
the question of their equivalence becomes non-trivial. In 
case of the standard BRM we have witnessed [|o| in the 
non-perturbative regime a premature departure from bal- 
listic behavior, and appearance of an intermediate dif- 
fusive stage. We observe essentially the same behavior 
in case of the EBRM Hamiltonian (Fig. Ic). But with 
the 2DW Hamiltonian (Fig. lb) we do not have such an 
effect: As ?i ^ the correspondence with the classi- 
cal behavior becomes better and better. Thus our sim- 
ulations demonstrate that having diffusion in the non- 
perturbative regime is an artifact of the RMT approach. 
In our previous work [|o| we did not have a numerical 
proof to support such a strong statement. There, all we 
were able to do, was to argue that RMT should fail in 
the deep semiclassical regime {fi ^ Csc), thus leaving 
open the possibility for having an intermediate regime 
(Csc < h < Cprt) where RMT might be valid. As we 



3 



-5 
-15 
-25 



-4 
-6 
-8 
-10 



2DW 

EBRM 

PRT 





-400 



-200 







200 



400 



FIG. 4. The spreading profile Poo{r) in representative 
cases. The upper panel is an example for a standard pertur- 
bative profile {P{t) ~ 1). The middle panel is an example for 
a perturbative core-tail structure {AxN <^ 5E). The lower 
panel is an example for an ergodic-like non-perturbative struc- 
ture (AxiV ~ 5^;). 



turbative core-tail calculation (PRT for brevity). The 
saturation profile is given by the expression P^(n\m) = 
X)n' \ {'n'{xo)W {x))\^\{n' {x)\m{xQ))\^ . It can be regarded 
as the auto-convolution of PE{n\m) = \{n{x)\m{xo))\^ . 
Thus the average saturation profile (r) is approx- 
imately related to the average local density of states 
Pii{r). The latter has been analyzed in 0. In Fig. 3 
we have calculated the PRT of P^{r) via an auto- 
convolution of the PRT of Pjiir). In the extended pertur- 
bative regime the major features of the saturation profile 
are captures by the PRT. The differences are mainly in 
the far tails where higher order perturbation theory is 
essential. There are also differences in the small scale 
details, where the non-perturbative mixing is important. 
In the non-perturbative regime the saturation profile be- 
comes purely non-perturbative, and the PRT becomes 
useless. This is because there is no longer separation of 
energy scales, which is the working assumption of the 
core-tail theory. 

We thank Felix Izrailcv for suggesting to study the 
2DW model. 



see, our numerical results do not give any indication for 
the existence of such intermediate regime. The failure of 
RMT happens as soon as we enter the non-perturbative 
regime {h < C^.t). 

As we make h smaller there is no indication in Fig. 1 for 
entering a non-perturbative regime. In we have 

made the important distinction between detailed QCC 
and restricted QCC. The former pertains to the whole 
spreading profile, while the latter pertains only to the 
variance. Restricted QCC is a robust type of correspon- 
dence that does not require a 'very small h\ If we want to 
have an indication for the crossover to a non-perturbative 
behavior we should look on other measures, such as N{t) 
of Fig. 2. If the spreading were classical-like, it would im- 
ply that the spreading profile is characterized by a single 
energy scale. In such case we would expect that N{t) 
and y/M{t) would be the same (up to a numerical fac- 
tor). Indeed this is the case for the Ti < 0.04 runs of 
Fig. 2. However, this is definitely not the case in the per- 
turbative regimes, where we have a separation of energy 
scales N{t) <^ \jM{t). In the perturbative regimes M (t) 
is determined by the tails, and it is insensitive to the size 
of the 'core' region. The width N{t) constitutes a practi- 
cal estimate for the latter (see Fig. 3a). It is N{t) = 1 for 
a stondarrf perturbative profile, and 1 <C N{t) <^ ^jM{t) 
for a fully developed core-tail structure. An alternate 
way to identify a standard perturbative profile is via the 
survival probability P{t). Indeed for Ti < 0.14 we see in 
Fig.3b that we have P{t) ~ 1. 

The difference between the perturbative and the non- 
perturbative spreading profiles is further illustrated in 
Fig. 4. Here we have plotted representative average sat- 
uration profiles, along with a comparison with the per- 
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